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ABSTRACT 

Clusters of galaxies are known to be dynamically active systems, yet X-ray studies of the low redshift population 
exhibit tight scaling laws. In this work, we extend previous studies of this apparent paradox using numerical 
simulations of two extreme merger cases, one is a high Mach number (above 2.5) satellite merger similar to the 
"bullet cluster" and the other a merger of nearly equal mass progenitors. Creating X-ray images densely sampled 
in time, we construct Tx, Mgas, and Yx measures within R^qq and compare to the calibrations of Kravtsov et al. 
(2006). We find that these extreme merger cases respect the scaling relations, for both intrinsic measures and 
for measures derived from appropriately masked, synthetic Chandra X-ray images. The masking procedure plays 
a critical role in the X-ray temperature calculation while it is irrelevant in the X-ray gas mass derivation. Mis- 
centering up to 100 kpc does not influence the result. The observationally determined radius R500 might conduce 
to systematic shifts in Mgas, and Yx which increase the total mass scatter. 

Subject headings: cosmology: miscellaneous - methods: numerical - galaxies: clusters: general - X-ray: 
hydrodynamics - scaling relations. 



1. INTRODUCTION 

The comparison between the observed mass function over 
cosmic time with the theoretical expectation constrains vari- 
ous cosmological parameters such as the total mass density, 
the dark energy equation of state and the normalization of the 

l2009l|Mantz et al.|2009| ). Such 



power spectrum ( Vikhlinin et al. 



a comparison requires a mode 



for relating observed cluster 



quantities to the underlying masses of the halo that support 
them. 

The works cited above use the intracluster medium (ICM) of 
clusters as a proxy for mass. Estimates of the total cluster mass 
can be derived from a hydrostatic assumption. Recent com- 
putational work indicates that such estimates are likely to be 
biased 10-20% low, on averag e, with a slightly larger disper- 
sion, relative to the true values ( Rasia et al.|2004 2006 INagai 
eFaL][2007l [Piffaretti & Valdarnini||20081 |Jeltema et al.||2008t 



Lau et al.|20Q9i r 



The dynamical youth of clusters implies incomplete thermal- 
ization of the gas, and some pressure support remains in the 
form of ICM bulk motions and turbulence. This result was 
seen in the first three-dimensional, gas dynamical simulations 
of cluster formation (Evrard 1990), and it appears to be inde- 
pendent of the numerical algorithm (Lagrangian or Eulerian) 
and of the physics so far encoded in simulations. On the last 
point, some caution is in order because non-thermal effects, 
such as cosmic-ray pressure and magnetic fields, that may af- 
fect the degree of turbulence (Lagana et al. 2009) have not yet 
been extensively studied. 

For cosmological studies, the slope, zero-point and scatter 
in the mass-observable relation can be introduced as additional 
degrees of freedom that are solved for in the li kelihood analysis 
([Majumdar & Mohr|2004] [Lima & Hu 20051 |Cunha & Evrard 



|2009| ). Prior constraints on these parameters, from observations 
or from simulations, can strongly enhance the statistical power 
of survey analysis. 

Scaling relations of the cluster population are thus an impor- 
tant tool for cluster cosmology. Recent computational modeling 
emphasizes the regularity of clusters, with the virial theorem a 
basic organizing mechanism. The calibration of dark matter 
(DM) velocity dispersion, (Jdm, by |Evrard et al.| ([2008) indi- 
cates that (Jdm at fixed mass and epoch is close to log-normally 



distributed with a scatter of4.3=b0.2%,a finding confirmed by 
the subsequent study of |Lau et"aL] ( |2009] ). From the observa- 
tional side, limits to using velocity dispersion as a mass proxy 
are imposed by projection effects, which substantially enlarge 
the scatter between mass and line-of-sight galaxy velocity dis- 
persion ( [Biviano et al. 2006), even when the latter is assumed 
to trace ctdm- Samples of roughly 4000 halos covering a range 
of redshifts from the Millenium Gas Simulations ( Stanek et al.l 



2009 [ [Short et al.|2010| ) support the log-normal assumption for 
X-ray observable signals, such as ICM mass Mgas, temperature 



Tx, and their product Yx ( [Kravtsov et al. 2006). The effective 
scatter in mass for these measures, while model-dependent, is 
15% or less for the case of a preheated ICM (Stanek et al.|2009| ). 

While such large samples naturally include systems that span 
a range of dynamical states, from dynamically old, relaxed ob- 
jects to dynamically young, merging systems, other studies have 
explicitly attempted to explore the role of mergers on scaling 
relations. [Ricker & Sarazin| (|20011) demonstrate the potential 
for large excursions in X-ray luminosity and temperature in 
the case of binary merg ers. This result contrasts with subse- 
quent work by Ritchie & Thomas] (|2002), who analyze differ- 
ent mergers scenarios including various impact parameters and 
two different mass ratios. These authors show that, while merg- 



ers do boost the X-ray luminosity and temperature (c.f. |Torri 
et aL\2004 ), they do not substantially increase the Lx-T scat 



ter. I Rowley et al.] ( [2004 ) confirm this result, finding no clear 
relationship between deviation about the mean Lx-T relation 
and merger activity or formation time. They report some evi- 



dence for a dependence on the degree of substructure. |0'Hara 



|et al.| (2006), using both simulations and observations, find no 
influence on scaling relations from morphological indicators of 



dynamical state. Finally, simulation analysis by Hartley et al 



( 2008| ) indicates that formation epoch is correlated with loca- 



tion in the Lx-T plane. They see weak negative curvature in 
the relation at the bright end caused by a prevalence of recently- 
formed systems. 

On the observational side, analysis of luminosity and en- 
tropy scaling relations in the REXCESS sample show that mor- 
phologically disturbed systems are offset relative to the whole 
sample, toward lower luminosities and higher entropies at fixed 
temperature (Pratt et al.||2009|). This appear s to be due to a 
deficit of gas interior to R500 ( Pratt et al.|2010| ). 
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While the scaling between X-ray luminosity and tempera- 
ture has a long history o f investigation (e.g.|Markevitch|1998] 



Amaud & Evrard|[T999l [Osmond & Po nman 2004; Prat t et al.| 
2009K the large scatter in the (non-core-excised) relation has 
led to the consideration of other mass proxies, especially the 
gas thermal energy, Yx, the product of the gas mass, Mgas, and 
temperature, Tx, derived from X-ray observations. 

Kravtso v"et al.| ( |20 06^ present the Yx parameter as mass proxy 
because of its low scatter and nearly self- similar behavior. Us- 



We summarize here only the basic concepts behind the sim- 
ulation process; further details can be found in |Dolag et"aL 



ing a small set of idealized binary merger simulations, [Poole 
et al.|(|200"7l) confirm the tight scatter of Yx as a mass proxy, and 



demonstrate that it is relatively insensitive to merger events and 
projection. With halos drawn from a large volume simulation, 
Yang et aL| ( |200 9 ) find that the scatter of mass-temperature re- 
lation depends on the halo concentration: at fixed mass, highly 
concentrated halos tend to be colder. That study did not find 
any bias between merging and non-merging systems. 

The simulation evidence cited above indicates that dynamical 
state plays less of a role than formation epoch in determining 
the scaling relation location of a particular cluster. Still, there 
remains concern that merging systems may bias cluster cosmol- 
ogy studies based on counts as a function of X-ray luminos- 
ity function or temperat ure ([Randall et al. 2002) or Sunyaev- 
Zel'dovich decrement ( Wik et al."2008V Here, we revisit the 
issue of the effect of mergers, emphasizing both theoretical 
and observational perspectives and using improved simulations 
with better resolution and more complex physics than used pre- 
viously. 

In this paper, we test the robustness of scaling relations in 
two specific merging systems. From a sample of 25 simu- 
lated halos ( [Dolag et al.||2009| ), we select two as examples of 
extreme cases which can occur in the universe. The first is 
a merger happening at z = 0.4 between a small object and a 
massive halo, M500 ~ 2 x lO^^M©. While the secondary struc- 
ture has a mass of only one-tenth this mass, its close passage 
through the core plasma causes a detachment of its baryonic and 
non-baryonic components. The Mach number of this powerful 
encounter is 2.5. Recent observations indicate that these rare 
circumstances are possible; witness the famous bullet cluster, 
1E0657-56, ( Ma rkevitch et al. 2002 ; Clowe et al. 2006) and the 
massive merger of MACS J0025.4-1222 ( [Bradac et aLl[2008l ) 
. We refer to this simulated system as the 'bullet-like cluster' 
(BL). The second situation is a merger happening at a more 
recent redshift (z ~ 0.09) between two halos of similar mass, 
^500 ~ 1.5 X IO^'^Mq. We refer to it as the '1:1 merger' case 
(1:1). 

With respect to previous studies, improvements in this work 
include the following: (i) our simulated systems evolve in a 
cosmological environment, therefore; (ii) our clusters experi- 
ence multiple interactions with the surrounding structures; (iii) 
the merging phases are followed with high temporal detail, and; 
(iv) we employ a rigorous X-ray spectral and imaging analysis 
instead of approximate formulae to calculate X-ray temperature 
and luminosity. 

The paper is structured as follow. Section 2 offers a descrip- 
tion of the simulated clusters. Section 3 describes the evolution 
of intrinsic cluster properties for both mergers. Results based 
on X-ray analysis are presented in Section 4, and sources of er- 
ror are dissected in Section 5. Section 6 offers our conclusions. 

2. SIMULATED CLUSTERS: SAMPLE 



([2009). Both clusters are originally identified in the same dark- 
matter-only cosmolo gical simulation, 479 Mpch~^ on a side 
( [Yoshida et al.[[2001| ), and then re- simulated using a more de- 
tailed treatment of the intracluster medium physics. The re- 
simulation technique used takes advantage of the Zoomed ini- 
tial condition method (Tormen et al. fl997) which largely in- 
crease the spatial and mass resolution. As result, the ' 1 : 1 merger' 
cluster has 2 million particles inside R200 Q while the 'bullet- 
like' system has a factor of 20 more particles. 

The re-simulations employ GADGET-2 ( [Springel[[2005[ ), in 
which we introduced a uniform and evolving UV background 
( [Haardt & Madau|1996[ ), star formation from a multiphase in- 
terstellar mediu m, a prescription for galactic winds triggered by 
SN explosions ("Spr ingel & Hemquist[[2003 ') and the influence 
of 1/3-Spitzer thermal conduction. In the 'bullet-like' simula- 



tion, we use the method by Dolag et al. (2005 ) to limit the arti- 
ficial viscosity typical of SPH simulations as well as the model 
of chemical enrichment by [Tomatore et al. p007| ). The cos- 
mology is a flat AC DM model with ^7o,m = 0.3 for the present 
matter density parameter, and ft^ = 0.039 for the baryonic den- 
sity, leading to a baryon fraction of 0.13. The Hubble constant 
is h = 0.7 in units of 100 km s~^ Mpc~\ and the rms level of 
density fluctuations within a top-hat sphere of 8 Mpc radius 
is (78 = 0.9. 

Bullet-like cluster. This halo is the most massive system of 
the 25 investigated clusters of the original sample (g^.a in Ta- 

^ R200 and R500 are defined as the radius of the sphere whose density is 200 and 
500 times (respectively) the critical density of the universe at that redshift. 
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Fig. 1. — DM (left panels) and ICM (right panels) particles of 
the substructure are shown in false color. The different rows corre- 
spond to the time tc (central panel), 600 Myr before (upper panel) and 
600 Myr after (lower panel). Two movies can be downloaded from 
[http://www.astro.lsa.umich.edu/~rasia/REL_SCA_MERGER.php 
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ble 2 of Dolag et al. 2009). It is in the central position of a 
very dense environment which includes four other systems of 
mass M200 > 1 X IO^'^Mq within a 10 Mpc radius. The halo is 
connected with its companions through several filaments from 
which it is accreting material. 

In Fig. [T] we show the dynamics of both dark matter (left 
column) and gas (right column) components of the substructure 
before, during, and after the merger. The substructure moves 
from North- West to South-East. In the top panels we show the 
moment (the fore ~ 8.5 Gyr, z = 0.49) when the two radii, R500 
(white circles), of the halo and of the secondary object intersect 
for the first time. The central panels show the instant (^^ ~ 9.1 
Gyr, z = 0.4) when the centers of the two systems are coincident. 
Finally, the lower panels show the situation (t after ~ 9.7 Gyr, 
z = 0.34) when most of the mass of the substructure exits from 
the 7?5oo sphere of the main object. The temporal range from 
tbefore ^nd tc and from tc and tafter is 600 Myr. 

During the interaction, the two components are experienc- 
ing different forces and are influenced in a different way by the 
merger. The gas is slowed down during the infall and shows 
already an asymmetrical tail at time the fore- During the merger, 
ram-pressure stripping confines the plasma to a region close to 
the center of the main cluster, and some gas stays even after 
tafter- Sincc the dark matter feels only the gravitational forces, 
after the passage through the center it continues its track rela- 
tively unperturbed. The star particle clumps that represent the 
galaxies in the simulation (not shown in the Figures) follow 
very closely the dark matter particles behavior, lagging only 
slightly behind the dark matter at time taper - 

While the most evident feature of this figure is the separation 
of the baryonic plasma from the non-baryonic component, other 
elements are distinguishable. First, the geometry of the cluster 
is not spherically symmetric, but elongated in West-East direc- 
tion. Second, note the presence of several other clumps which 
fall into the main halo. Among them a chief role is played by 
a very cold and massive satellite which is located on the West 
side on the central panels and on the South- West on the lower 
images. In what follows, we denote as tcoid the instant when 
this satellite enters within Rsoo of the main cluster. 

To trace the interaction evolution we use 26 snapshots from 
redshift ~ 0.5 (or Hubble time equal to ^ = 8.2 Gyr) to redshift 
0.28 (t = 10.2 Gyr). The time step between consecutive snap- 
shots is generally 150 Myr, excluding the period from tc and 
tafter whcrc a shorter interval of 30 Myr is used to increase the 
temporal resolution. 

1:1 merger case. The second situation we follow is the final 
phase of the formation of a moderate mass halo with M200 = 
5 X 10 ^^M0 at z = 0. As in the former case, this halo also lies 
in a dense environment including 5 objects with M200 greater 
than 0.5 x 10^"^ within 10 Mpc. However, it is not the most 
massive member of this supercluster (gl.b in Table 2 of Dolag 
et al. 2009) 

We study the most recent major merger involving two objects 
of comparable mass. In Fig|2j we show 16 intrinsic flux maps 
that exhibit the evolution of this merger. The lower mass sys- 
tem arrives from South-East and passes to the East side of the 
center. It reaches a closest distance of 300 kpc at time t^ and 
executes a spiral motion before being complete absorbed by the 
principal halo. The line of sight that we chose is that which 
minimizes the impact parameter and might create the biggest 
mis-interpretation of the X-ray data. Indeed, the flux map cor- 




FiG. 2. — Intrinsic flux images of the '1:1 merger case' taken from redshift 
0.3 to redshift 0.06. The field of view is fixed to 4 Mpc. The color scale is the 
same for all panels and is in units of photons/s/cm^. 



respondent to z = 0. 1 1 in Fig. 2 present regular iso-flux contours, 
indication of a relaxed cluster, and only one bright spot. 

For the rest of the paper, we consider 15 snapshots from red- 
shift 0.3 (t = 10.07 Gyr) to redshift 0.09 (t = 12.23Gyr), sepa- 
rated by 150 Myr. As above, the fore denotes the instant when 
the two Rsoo are in contact. 

3. INTRINSIC SCALING EVOLUTION 

For each snapshot, we locate the center of the halo as the 
minimum of the potential well. We compute the gas mass, 
M gas, Sim, Summing the gas particles inside R500, and the total 
mass, Mtot, considering the contribution by all species of parti- 
cles (dark matter, gas, and stars) . The choice of this radius is 
motivated by both t heoretical ( [Evrard et al.|1996[ [Rowley et al 
20041 Nagai et a L|r2007) and observational works ( Vikhlinin 
et al.||2006 ) that establish this as a scale within with scaling 
relation scatter is minimized. 



Table 1 

Parameters of the fitting formula 

Mu^tot'^Eiz)^ = Mo(X/jCo)" IN THE CASE OF ALL REDSHIFTS AND 

ALL CLUSTERS (TAB. 2 0F |KRAVTS0V ET AL7r 2006). Mq IS 
CHANGED ACCORDINGLY TO THE DIFFERENT BARYON FRACTION 



Mo 



Xo 



scatter 



r[keV] 

MgasllO'^M^] 
YxIIO^'^MqX keV] 



2.55 
2.42 
1.95 



1.521 
0.921 
0.581 
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0.4 



19.5 
10.7 
7.1 
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To identify the temperature, we apply the spectroscopic-Hke 
definition, 7^/ = ^WT dV / dV with T the temperat ure of 
each gas pa rticle, n its density and W equal to n^/T^^-^^'' (Maz-| 
zotta et aL|2 QQ4 ). The sum is extended to all the particles with 
temperature larger than 0.5 keV within [0.15 1] R^oq. The ex- 
clusion of the inner region is conventional among observers and 
simulators since it guarantees that X-ray results are not influ- 
enced by the activity of the cluster core and intrinsic analysis 
of simulations are not affected by the overcooling problem. Fi- 
nally, Yx^sim is the product of the previously measured quantities 

In Fig. [3] we show the time evolution of the total mass, the 
Yx.sim parameter, the gas mass, and the spectroscopic-like tem- 
perature of both halos. The evolutionary tracks of the masses 
of the bullet-like halo show several bumps which testify to the 
continuous accretion of material due to the infall of satellite ob- 
jects (see also, Fig.[T]). Notice that many of these show a mass 
ratio which is higher than that we are focusing on in this work 
(1:10). Nevertheless, this particular merger is the only one that 
produces an increase on temperature of 40% and a Mach num- 
ber of 2.5. After tcoid^ the temperature has a rapid drop due to 
the arrival of the second, colder satellite. The decline is shal- 
lower for Yx.sim due to the combined effect of the increase of 

^gas^sinf 

In the case of the '1:1 merger' the evolutionary tracks are 
smoother and they all present only one big bump of order of 
100% for the masses and 40% for the temperature due to the 
merger. In both cases, the increase in the temperature is not 
simultaneous with the increase in the mass, a delay of few hun- 
dred mega years is needed to allow shocks to develop and heat 
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Fig. 3. — Evolution of the spectroscopic-like temperature, Tsi[keV\ gas 
mass, Mgas^simVl^^^MQ'], Yx^sim =Mgas,sim X Tsi[1()^^Mq X keV] parameter, 
and total mass, M^o^[lO^^M0] (bottom to top panels). On the left side we 
present the 'bullet-like' case, while on the right side the '1:1 merger' case. The 
vertical lines indicate important moments in the clusters' evolution: f^g/ore, tc, 
and tcoid (see text in Sec|2]for details) 



the intra-cluster medium. 

4. X-RAY-OBSERVED SCALING EVOLUTION 

We create X-ray images for each snapshot with our X-ray 
MAp Simulator (X-MAS). The characteristics of this software 
package are described in detail in other works ( [Gardini et al. 
20041 [Rasia et al.||2006l [2008] ). To create the photon event 
file, we use the ARF and RMF of the Chandra ACIS-S3 de- 
tector aimpoint. We fix the redshift at the value of the simu- 
lated time frame and the equivalent hydrogen column density to 
Ufj = lO^^cm"^. In the '1:1 merger case', the metallicity is one 
third of the solar abundances by ,Anders & Grevesse[ ( [1989] ). In 
the 'bullet-like' halo, instead, we use the de tailed star forma- 
tion history model by [Tornatore et al.[p00 7) and each particle 
has different abundances for C, N, O, Mg, Si, Fe and a constant 
value for the remaining elements (see also Rasia et aLl[2008[ ). 
The field of view of our images is 16 arcmin and the depth we 
use in the projection is 10 Mpc. The exposure time is 1 Msec. 

Selected photon images in the soft band, [0.7 2] keV, are pre- 
sented in Fig. |4] for both clusters. In the images of the bullet- 
like cluster, the position of the bullet and of the cold object 
discussed in Sec. [2] are evident. As in Fig. [3 we note also the 
presence of many other clumps with comparable size to the bul- 
let. All these structures play an important role in perturbing and 
modifying the ICM appearance. The '1:1 merger' cluster also 
shows the existence of small blobs, but they are much less mas- 
sive than the merging object we are considering. 

For each time step, we apply the following X-ray spectral 
image processing to obtain the gas mass Mgas,x , the tempera- 
ture Tx and the parameter Yx. It is important to stress that we 
proceed on every single time step in an independent way, i.e. 
without knowledge either of the information stored in the sim- 
ulation itself or of its past history or future evolution. 



Surface Brightness. Using CIAO tool (Frus cione et al.[2006| ) 

we extract soft band images in the [0.7 2] keV band per each 
snapshot and use them to identify cold dense regions with the 
wavelet algorithm of [Vikhlinin et al.[ ( [l998[ ). These regions are 
masked and excluded from all the analysis. We proceed to ex- 
tract the surface brightness profiles in 30-40 linearly spaced an- 
nuli. For the center, we generally use the centroid of the exter- 
nal and more spherical iso-flux contours. The profiles, stored as 
counts versus radius, span over a radial range from ~ 0. 1 x R^qq 
to - 2 X R500. 

Temperature. To measure the X-ray temperature, we extract 
the spectrum in the [0.5-7] keV band of the region inside R500 
with the exclusion of the masked regions and of the inner 15% 
of R500. At a first pass, R500 is defined directly from the simu- 
lation (see Sec|2]). In a second step, we derive R500 following 
the X-ray procedure described below. We refer to this measure 
as Rsoo^x- The X-ray spectra are fitted using a statistic in 



the XSPEC package ( [Arnaud[[1996[ ) with a single temperature 
MEKAL model by fixing input values of the hydrogen column 
density (hh = lO^^cm"^), redshift and, in the '1:1 merger case', 
the metallicity (Z = 03Zq). 



Gas Mass. Following Vikhlinin et al. ( [2006 )'s procedure, we 
fit simultaneously the surface brightness and the temperature 
profiles using the following expressions: 
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Fig. 4. — Photon images in the [0.7 2] keV soft band of 'bullet-like' (left panels) and '1:1 merger' cases (right panels). The field of view is 16 arcmin and the 
exposure time 1 Msec for each panel. The entire collection of X-ray images can be found at http: //www.astro.lsa.umich.edu/^rasia/ANALYSI S_bullet.php and 
http://www.astro.lsa.umich.edu/~rasia/ANALYSIS_l_l .php 
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The surface brightness fitting formula is an extension of the 
P model ( [Cavaliere & Fusco-Femiano|1976| ) with the addition 
of an inner power-law cusp (the a term), and a change in the 
slope in the external region (slope modified by e around r^). The 
temperature expression is the product of a broken power law. A 
more comprehensive description includes additional terms that 
describe the inner core region (see |Vikhlinin et al.|[2QQ6| ), but 
we do not consider them here since we exclude the inner 15% 
of Rsoo from our profile analysis. 



■500,X« 



We derive an estimate of the scale radius, R 



500,X, US- 



ing only the X-ray information. This iterative method ( [Kravtsov 
et al.||2 006 ) is based on knowledge of the M-Yx scaling rela- 



tion, of the gas mass profile, Mgas, and of different measure- 
ments of the temperature, Tx, inside subsequent estimate of 
R500,x- Calling the first guess OLD, the equation we need to 
solve to find the NEW Rsoo^x is based on the two definitions of 

Mtnt. 



47r 



^loO^X^NEWPcr = MoE(zf (Yx,OLD / XqT 



= MoEizfrn^as^OLD X rx,OLD) Aol", (3) 

where pcr is the critical density at the considered redshift, z; 
E(z) = H(z) /Ho = [(^m(1 +1^a]^-^ (Peebles 1993") is the evo- 
lution factor; Mq, xq, a and /3 are the specific constants of the 
M-Yx relation (see Table[T]). Yx^old, Mgas.OLD, and Tx^old are 
the quantities measured inside the previous attempt of Rsoo^old- 
This equation defines the next value, Rsoo^x^new^ for which 
we repeat the measurements of temperature and the calculation 
of gas mass. The process is iteratively repeated until subsequent 
estimates of R^oo are consistent within 3%. In our sample, we 
begin our calculation assuming the radius computed from the 
simulation directly, Rsoo.sim- For this reason, in 75% of the cases 
we stop the iteration after the first round. 
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Fig. 5. — Scaling relations from intrinsic analysis of the simulated halo. Left 
panels refer to the bullet-like halo, while right panels to the 1:1 merger case. 
From the top, we show Mtot -^gas.simi^tot -Tsi,andMtot -Yx^sim- The total 
mass is multiply by E(z) to a power equal to 0, 1, 2/5, respectively (see text 
for the definition of E(z)). The different color are the same ofFigp^ Red solid 
and dashed lines represent Kravtsov et al. (2006) relation and its scatter (see 
Table.|2). 
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5. RESULTS 

The scaling relations results for the quantities intrinsically 
computed in the simulation (SecjS]) and the measurements done 
following the X-ray procedure (Sec. |4]) are presented in Fig. |5] 
and Fig.[6j respectively. In both figures, the j-axis is to the total 
mass multiplied by the evolution factor, E(z) at the appropriate 
f3 power (see Table 1). Evolution powers are a = for gas mass, 
a = 1 for temperature, and a = 2/5 fovYx. The color scheme for 
points in Figs. [5] and |6] is the same used in Fig[3] with the black 
symbols specifying the significant moments, t before^ tc^ and, in 
the case of the 'bullet-like' cluster, tcoid- 



In every panel, the solid red line shows the Kravtsov et al 



( |2QQ6| ) scaling relation rescaled by the mean baryon fractions 



of our respective simulations. We assume a baryon fraction 
VLt = 0.039 (see Sec|2]) while they assume = 0.04286. The 
ratio, 1.099, corrects Mgas- Considering that their average < 
Mgas/Mtot > was ~ 0.09, we correct also Mtot by 1.008. Fi- 
nally, the dashed lines represent their estimate of the la scatter 
in each scaling relation. Tab[T] lists parameters of the scaling 
relations shown, rescaled for our cosmology. 

5.1. Intrinsic scaling relations 

The two systems under study have substantially different evo- 
lutionary behavior; the 'bullet-like' cluster has a complicated 
history with the multiple accretion events while the '1:1 cluster' 
evolves more smoothly. The BL variation of mass and tempera- 
ture are not synchronized, so there are several loops in the Mfot - 
Tsi relation. Moreover, the presence of the secondary, cold ac- 
cretion causes significant movement in the Mtot - Tsi plane to- 
ward the cool side of the mean relation (transition between the 
cyan points before tcoid and the green points after tcoid)- This 
variation is also reflected in the Mtot-yx,sim relation. 

The scaling evolution of the 1 : 1 case is more regular. Except 
for a phase delay during the merger itself, the gas mass, temper- 
ature, and Yx^sim observables (top to bottom) in Fig. [5] increase 
in a manner that largely reflects the appropriately scaled growth 
in total mass. The temperature and Yx^sim parameter experi- 
ence delays of ~ 150 Myr (see also Fig. [3]), as shocks expand 
and heat the plasma. That produces jumps in the Mtot - Tsi and 
^tot - Yx,sim relations that changes the sign of the mass proxy 
residuals, from slightly underestimating to overestimating the 
true mass. Yet the entire evolutionary histories of the 1:1 halo's 



scaling relati ons lie within the l-a region proposed by |Kravtsov 
et al.| ( 2QQ6| ). Since their analysis was based on measures de- 



rived from synthetic X-ray observations, it is not surprising that 
our intrinsic data should show smaller scatter. We turn now to 
our estimates of X-ray observable quantities. 

5.2. X-ray observed scaling relations 

The scaling relation between the total mass and the gas mass 
in Fig|6] shows a similar behavior of Fig|5j the two systems 
move along the [Kravtsov et al.| ( |200 6 ) relation without major 
deviation caused by the merging. In general, the X-ray gas 
mass estimates are slightly larger than the true gas ma ss, by 
5.5% on average. This value agrees with the results of | |Nagai| 
|et al.| ( |2QQ7 ), who found a bias of 5.8% ± 5% for their unrelaxed 
systems at redshift zero. The X-ray temperatures are hotter for 
the bullet halo in the final phase of the merger and, in particular, 
after tcoid (green points). The temperatures of the 1:1 case track 
the intrinsic evolution with a high degree of fidelity. Finally, the 
combined behavior of the gas mass and temperature lead to Yx 
being slightly higher in the X-ray analysis. 
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5.3. A priori difference between the two procedures. 

The differences between the intrinsic and X-ray scaling re- 
lations are slight, and some of the difference can be explained 
by bias in setting the region within the quantities are measured. 
In this section, we examine whether and how the results are 
affected by the choice of the center, the radius and the masking. 

Center. In the intrinsic analysis of the simulated halos we 
define the center as the minimum of the potential well while in 
the X-ray analysis the center coincides with the X-ray centroid. 
The two definitions agree well for the 1:1 case, but they often 
differ in the BL case, by typical values of 90 kpc. Such an off- 
set introduces a negligible effect in the temperature calculation. 
While it might change the inner profile of the surface bright- 
ness, and thus modify the value of the gas mass, our analysis 
masks the innermost regions where this effect would be rele- 
vant. To explore the importance of the center definition, we 
test whether there is a correlation between the centers' offset 
and the gas mass difference and found only a negligible posi- 
tive correlation ( Pearson correlation coefficient equal to 0.3). 
We conclude that the gas mass difference cannot be explained 
by the different selection of the center. 

Radius. On average, the true radius and the X-ray one are 
close. The maximum radii deviation. A/? = (^5oo-^500,x)/^500. 
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Fig. 7. — Am in colored points, /S.r is black diamonds and /S^Mxr in red. 
See text for the definition of the various A. 



Table 2 

Percentage deviation between the true radius and the 
x-ray radius, |ai?|, the intrinsic gas mass and the x-ray 
ONE, |Am|; the X-ray gas mass at the true radius and at 

THE X-RAY RADIUS, \/\mxA'^ THE X-RAY GAS MASS AND THE 
INTRINSIC ONE AT THE TRUE RADIUS, |AM,i?5oo|; AND THE 
CORRELATION BETWEEN /\r AND Am- 



BL 



1:1 



|A(7?)| 1.8% 3.7% 

|A(M)| 6.2% 8.1% 

|A(Mx,7?)| 1.7% 4.7% 

|A(M,7?5oo)| 4.9% 3.7% 

correlation(A/?,AM) 0.70 0.95 



is driven by the merger and it is always below 5%. However, 
the sign of the deviations differ: the X-ray radius is greater than 
the true one after the merger for the 'bullet- like halo' , while 
it is the opposite for the '1:1 merger case'. We remind that in 
the 7?5oo,x computation, the 3% requisite for further iterations 
is applied to two subsequent radius evaluations (see Sec. 4), 
the final difference with the intrinsic radius might, therefore, be 
higher. 

The small radial discrepancy directly affects the X-ray gas 
mass estimate since, at 7?5oo, the gas density profile has a slope 
close to -2, implying that the gas mass is proportional to R 
( Vikhlinin et al.|20()6| ). A deviation of few percent in the radius 
implies a similar fractional error in gas mass. 

In the bottom left panel of Fig|7] we sketch the situation. The 
two dashed-dotted lines represent the intrinsic (lower) and es- 
timated X-ray (upper) gas mass profiles. The two large points 
indicate the value of the two gas mass at the two respective 
radii. The difference between the two radii, A(/?) (the two hor- 
izontal lines), is similar to the difference of the X-ray gas mass 
computed at the two radii, A(Mx,7?) (short red vertical line at 
the right). For our halos, the situation is outlined in Table [2] and 
on the right panels of Fig|7j the A (7?) shifts (black diamonds) 
almost always coincide with the A(Mx,7?) (red crosses). The 
total mass difference, A(M) (the long blue vertical line at the 
right in the bottom left panel) is represented by colored points 
in the left panels. In addition to A(Mx,7?) it includes also the 
discrepancy between the X-ray and the intrinsic gas mass com- 
puted at a fixed radius (dashed vertical line on the bottom left, 
A(M,7?5oo)). The last difference could be caused either by too 
restrictive assumptions on the geometry of the systems or by 
the presence/exclusion of the cold blobs. Even if these effects 
are hardly controllable, we can still reduce the systematics by 
computing correctly the radius. This will decrease the disagree- 
ment between the two masses from 6.2% to 4.9% for the BL and 
from 8.1% to 4.7% for 1:1 case. Notice that the X-ray radius 
is both an under- and over-estimate of the true value. However, 
the implication on the gas mass are larger in the over-estimate 
scenario. As consequence a wrong X-ray radius derivation does 
not simply increase the scatter but skews the gas mass distribu- 
tion, changing, de facto, the scaling relation. 

These variations in mass are consistent with the results of 
Meneghetti et al.|([2009) . Using simulations from the our sam- 
ple ( Dolag et al.|2009 ) they found a difference in gas mass of 
7% ± 3% when computed at the same lensing radius. 

The Fx -based scaling technique to compute the mass and ra- 
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dius, described in Sec.|4] has inherent scatter that drives system- 
atic error. If the cluster is naturally below the Mtot - Yx scaling 
relation (such as the 1 : 1 halo at the beginning), the X-ray radius 
inferred through Eq. 3 will be greater than the true one. As con- 
sequence the cluster size will be shifted towards larger values 
in both the Mgas and Yx measures. The increase of the scatter 
along these proxies-axes will be the same, however, the total 
mass scatter will be smaller for Yx which has a lower power 
law. 

Masked regions Another source of difference from the pro- 
cedures followed in Sec.|3]and Sec.[4]is the masking technique. 
In the X-ray analysis we use a wavelet algorithm to exclude lo- 
calized peaks present in the intra-cluster medium, and to com- 
pute the spectroscopic-like temperature we cut all the particles 
below 0.5 keV independently on the particle density or position. 
Masking the cold regions is a fundamental step to link our work 
to the observed X-ray clusters. With this procedure we ensure 
that (i) the results based on the mock images are not influenced 
by the presence of the overcoolin g problem affecting radiative 
simulation (Bo rgani & Kravtsov||2 QQ9 ); (ii) the appearance of 
the remaining X-ray emission is more similar to real X-ray ob- 
servations (Boh ringer et al.||2Q10| ); (iii) the X-ray analysis we 
perform mirrors observers approach. 

The measurement most affected by the masking is the tem- 
perature, as comparison of Fig|5]and Fig |6] attests. In the X-ray 
analysis, we exclude the secondary cold satellite. As a con- 
sequence, the X-ray temperature between tc and tcoid (green 
points) does not drop as dramatically as the intrinsic spectroscopic- 
like one, and the mass -temperature scatter is therefore lower. 

6. CONCLUSION 

We study the behavior of X-ray proxies for total mass (gas 
mass, temperature, and their product, Yx) in case of two extreme 
merger events. The evolution of intrinsic, three-dimensional 
measures is compared to estimates derived from projected, X- 
ray spectral maps with high signal-to-noise. We create photon 
event files using our X-ray MAp Simulator and analyze them 
with standard procedures. 

Specific behaviors of the scaling relations, and lessons learned 
from comparing intrinsic measures with those determined from 
synthetic observations, are: 



• Code comparison. After a small correction for cosmic 
baryon fraction differences, the scaling relation behav- 
ior of both merger scenarios lies within the ±2<j regions 
expected from ART simulations ( [Kravtsov et al.|2QQ6| ). 
This agreement in non-trivial, given our use of a differ- 
ent code (GADGET) employing a somewhat different 
physical treatment. 



• Gas mass. The intrinsic gas mass is a very good proxy 
for the total mass, growing almost contemporaneously 
and linearly with it in both merger cases. The X-ray 
gas mass estimate is biased high by typically 5%, due to 
asphericity in the hot ICM that is not accounted for in 
the model. The gas mass measurement is not affected by 
the masking technique but strongly depends on possible 
bias in the radius definition (see below). 

• Temperature. The temperature is the X-ray proxy that 
shows the largest scatter and is the most influenced by 
the presence of substructures in the intra-cluster medium. 



Masking cold satellite emission reduces the scatter sub- 
stantially. 

• X-ray SZ Estimate, Yx» The scatter in the mass-l^ re- 
lation is maintained low during the both entire merging 
processes. 

• Center. The cluster center is defined independently in 
the intrinsic and X-ray analysis, being the minimum of 
the potential well in the former and the X-ray centroid in 
the latter. Mis-centering differences up to ~ 100 kpc are 
present for the BL cluster, but they do not significantly 
influence the mass proxy estimates at large radius. 

• Radius. Estimates of the radial scale, 7?5oo, do affect 
mass proxy estimates, with typical error 5% in both 
merger cases . The X-ray iterative derivation of the ra- 
dius might conduce to an over-estimate with reflections 
on the gas mass measurements. The total mass scatter 
of the Mtot -Mgas scaling relation will be more affected 
than that of Mfot-Yx due to the different power laws. 

• Masking. The masking procedure plays a critical role 
when comparing radiative simulations with X-ray ob- 
servations. Results from simulation that do not exclude 
any regions can overestimate temperature deviation with 
respect to the mean. We show that this has a large im- 
pact on the temperature calculation but not on the gas 
mass estimate. 



The results of this paper, combined with the finding by |Yang| 
et al.| ([2009) that the scatter is lower at higher redshift even con- 
sidering merging clusters, are encouraging for future missions 
(such as e-ROSIT.^or the Wide-Field X-ray Telescop^ that 
aim to observe clusters at high redshift, where more frequent 
mergers are expected. 
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